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Abstract: We report a mathematically rigorous technique which facilitates 
the optimization of various optical properties of electromagnetic fields in 
free space and including scattering interactions. The technique exploits the 
linearity of electromagnetic fields along with the quadratic nature of the 
intensity to define specific Optical Eigenmodes (OEi) that are pertinent to 
the interaction considered. Key applications include the optimization of the 
size of a focused spot, the transmission through sub-wavelength apertures, 
and of the optical force acting on microparticles. We verify experimentally 
the OEi approach by minimising the size of a focused optical field using a 
superposition of Bessel beams. 
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1. Introduction 

The decomposition of fields into eigenmodes is a well established technique to solve various 
problems within physical sciences. The most prominent example is the Schrodinger's equa- 
tion within the field of quantum mechanics, where energy spectra of atoms are determined via 
the eigenvalue spectra and associated wavefunctions of the Hamiltonian operator. Indeed, elec- 
tron orbits are eigenmodes of the energy, angular momentum, and spin operators (I) and as 
such they deliver fundamental insights into the physics of atoms. Within classical mechanics, 
modes of vibration of music instruments give, for example, their resonant frequencies while 
their spectrum is associated with the shape of the instrument [2]. In the optical domain, mode 
decomposition is used in order to describe light propagation within waveguides [3 ], photonic 
crystals [4], optical cavities 0, laser resonators 0, and the optical forces on Mie-sized par- 
ticles [7 ]. In the case of waveguides and photonic crystals, for example, eigenmodes describe 
electromagnetic fields that are invariant in their intensity profile as they propagate along the 
fibre or crystal. Additionally, these modes are orthogonal and as such light coupled to one of 
these modes remains, in theory, in this mode forever. This optical mode decomposition can be 
expanded to include additional properties such as orbital and spin angular momentum 0. All- 
together, the eigenmode expansion method is a well-established method for the representation 
of the propagation of optical fields. 

In this paper, we report a novel method which we term "Optical Eigenmodes (OEi)" which 
represents a generalization of the powerful concept of eigenmode decomposition going beyond 
the propagation properties of light. Crucially, we show that eigenmode decomposition is appli- 
cable to the case of any quadratic measure which is defined as a function of the electromagnetic 
field. Prominent examples of optical quadratic measures include the energy density and the en- 
ergy flux of electromagnetic fields. The OEi method makes it possible to describe an optical 
system and its response to incident electromagnetic fields as a simple mode coupling prob- 
lem and to determine the optimal "excitation" for the given measure considered. Intuitively, a 
superposition of initial fields is optimized in a manner that the minimum/maximum measure 
is achieved. For instance, the transmission through a pinhole is optimized by maximizing the 
energy flux through the pinhole O. 

From a theoretical perspective, the OEi optimization method is mathematically rigorous and 
may be distinguished from the multiple techniques currently employed ranging from genetic 
algorithms [10] and random search methods [11] to direct search methods fl2l . The major 
challenge encountered in any such approximate optimization and engineering of optical prop- 
erties is the fact that electromagnetic waves interfere. As such the interference pattern not only 
makes the search for an optimum beam problematic but crucially renders the superposition 
found unreliable, as the different algorithms may converge on different local minima which 



are unstable with respect to the different initial parameters of the problem. In contrast, our pro- 
posed OEi method yields a unique solution to the problem and directly determines the optimum 
(maximal/minimal) measure possible. 

In the first part of the paper, we introduce the OEi method and show its properties in a 
general context of optimizing the quadratic measures of interfering waves. In the second part, 
we apply the OEi formalism to minimize the focal spot size and discuss the appearance of 
superoscillating fields. For these applications, we describe, respectively, the electromagnetic 
field as a superposition of scalar Laguerre- Gaussian beams, vectorial Bessel beams or more 
general plane waves within the angular spectral decomposition representation of light. In the 
third part of the paper, we report a particular experimental implementation of the OEi method 
using computer controlled spatial light modulators to squeeze the spot size of a superposition of 
Bessel beams. In the last part of the paper, the method is applied in numerical 3D modelling to 
determine the OEi yielding the largest transmission through a sub- wavelength aperture and the 
largest optical force on a micrometer sized particle. The paper concludes with a discussion of 
the particular results obtained and with general comments on the versatility of the OEi method 
to a wide range of problems. A short annex compares the convergence properties of the OEi 
approach with standard phase front correction methods. 

2. Method 

The OEi method exploits both the linearity of Maxwell's equations and the quadratic depen- 
dence of light-matter interactions on the electromagnetic field {E, H} where E and H denote the 
electric and magnetic field vectors, respectively. Table [T] provides a list of common examples 
of such interactions. These interactions may be written in a general quadratic matrix form 

mW(E,H)=a t M( A 'a (1) 

where we considered a superposition of fields {E,H} = ^Y!J=i a j iY!J=\ a j^- an d where 

(A) labels the light-matter interaction defined in Table [l] The vectors a and are comprised 
of the superposition coefficients aj and their complex conjugates, respectively. The elements 

Mj^ of M( A ) are constructed by combining the respective fields {E^H,-} and {E^,H^} for 
j,k= 1 . . . N. More precisely, we have: 

4M^ } = m( A )(E i -+E^H^ + H^)-/m( A )(E i - + /E^H^ + /H^) 

-m^ (Ej - E h H k - H*) + irnW (Ej - iE h H k - iU k ) . (2) 

Given the Hermitian form of ([2]), we remark that the light-matter interaction M^ A ) defines 
a spectrum of real eigenvalues xjf^ and associated eigenvectors v[ A ^. Each of these eigen- 
vectors corresponds to a superposition of fields {E^Hy} termed here Optical Eigenmode 
(OEi). Crucially, we may now extremize the light-matter interaction considered; that is, the ex- 
tremal eigenvalue X^} and the associated eigenvectors v^t deliver the superposition of fields 
{E ext , H ext } = |l7=i v ext > ^j=i v ext j^j } wrncn extremizes the interaction while keep- 
ing the total field contributions constant a^a = 1. 

In this paper, we apply the OEi concept to minimize the size of a laser spot within a surface 
ROI. One way to define the spot size of a laser beam is by measuring (whilst keeping the total 
intensity constant) the second order moment w of its intensity distribution fT3ll . Crucially, w can 
be expressed in terms of and (see Table[l]) or the respective matrix representations ([]} 
as follows: 

„ /mP) „ /atM( 2 )a 



where and M^ 2 ) are termed the intensity operator (10) and spot size operator (SSO), 
respectively. According to ([3]), the minimum spot size is obtained by the OEi associated with 
the smallest eigenvalue A^ 2 J of the SSO provided that the 10 is simultaneously diagonalized 
and normalized to 1 . Direct evaluation shows that this is precisely achieved by the combined 
OEi 

f N N v ( 2 .) J?} N N v ( 2 .) JO) ) 

{BW,IW =11 -^7^ - E,, Z Z ^7="^ ■ < 4 > 

where vffi n k is the eigenvector associated with the smallest eigenvalue of M^ 2 ) and in the inten- 
sity normalised eigenbase v[°j / yj of M (°) . 



Energy 


mW(E,H) = 




Intensity 


m(°)(E,H) = 


i JT 5 (E* x H) • nda+c.c. 


Spot size 


m( 2 )(E,H) = 


\ J s r 2 (E* x H) • ndcr+c.c. 


Momentum 


m( F ' u )(E,H): 


= J / 5 (£b(E* • n)E + Mo(H* • n)H - \Sn) • u rfa+c.c. 



Table 1. Time averaged quadratic measures m of common light-matter interactions. The 
integration either over a volume V or a surface 5 which in general corresponds to the 
Range of interest = ROI of the measure. In the optical momentum case, it corresponds to 
a closed surface surrounding the scattering object with F • u representing the optical force 
in the direction defined by the unit vector u. For surface integrals, n is the normal unit 
vector to the surface considered. The definition of the electromagnetic energy density is 
$= l/2(e E-E*+ J u<)H-H*). 



For our proof-of-principle studies described in the remainder of the paper, we also applied 
the scalar version of the OEi method where a set of scalar fields E[ is considered in order to 
determine the 10 and SSO as 

Mf=jE]E k da (5) 

and 

M$ = fr 2 E]E k da, (6) 

respectively. These scalar expressions are equivalent to the respective vector versions listed in 
Table [T] and determined through equation ([2]). The scalar version of the optimimum OEi §5§ 
explicitly reads 

N N v ( 2 .) JO) 

E mn =^^^M-Ej. (7) 

j=lk=i ^ (0) 

2.7. Smallest focal spot using Laguerre Gaussian beams 

Using a superposition of LG beams we can minimize the size of a focal spot using the repre- 
sentation of the SSO ([6]) in cylindrical coordinates. It is important to note at this point that we 
only retain the intensity OEis whose eigenvalues are within a chosen fraction of total intensity. 
This is equivalent to considering only beams that have a significant intensity contribution in 
the ROI. Intuitively, the optimization procedure may perform so well that a spot of size zero is 
finally obtained if no intensity threshold is applied. Figure [T] shows the smallest spot superpo- 
sition where we observe the appearance of sidebands just outside the ROI. These sidebands are 
a secondary effect of squeezing the light below its diffraction limit. It is these sidebands that 



decrease the efficiency of the squeezed spot with respect to the maximal possible intensity in 
the ROI as calculated via the 10. Using the ratio between these two intensities we can define 
the intensity Strehl ratio ifTH for the SSO (see Fig. |2j)). We remark that both, the spot size 
and the Strehl ratio, show resonances as a function of the ROI size. This can be explained by 
considering the number of intensity eigenmodes used for the spot size operator. Indeed, as the 
ROI size decreases, so does the number of significant intensity eigenmodes. Each time one of 
these modes disappears (step in Fig. [2]), we have a sudden increase in the minimum spot size 
achievable accompanied with an enhanced Strehl ratio as we drop the most intensity inefficient 
mode. Overall, the Strehl ratios determined in our studies predominantly exceeded values of 
1% even when spots were tightly squeezed. Therefore, the observed decrease of intensity is not 
to severe in terms of potential applications of squeezed beams such as optical manipulation and 
imaging. 
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Fig. 1. (a) Transversal and (b) longitudinal 2D intensity cross sections of the OEi super- 
position delivering the smallest focal spot in the ROI (R = A) considering 25 LG modes. 
w/wq is the relative spot size measured according to Eq. ([3}. The Strehl ratio in (a) is 4.5%. 
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Fig. 2. (a) Spot size as a function of the radius of the ROI for different number of LG 
modes considered. The right hand scale and filled curve indicate the numbers of intensity 
eigenmodes fulfilling the intensity criteria for the N = 11 case. The arrows indicate 
the corresponding scales, (b) Ratio between the ROI intensity of the smallest spot size 
eigenmode and the largest intensity achievable in the ROI (Strehl ratio). 



On a final note, we remark that squeezing light below its diffraction limit may be associ- 
ated with the effect of super-oscillations fT5ll . This refers specifically to the ability to have a 





-15 -10 -5 5 10 15 
Wavevector (2n/wo) 



-I 



-2 




-2 



-1 





x/Wq 



Fig. 3. (a) Radial wavevector spectral density. Yellow highlights regions outside the spectral 
bandwidth, (b) Transversal cross section of the OEi spot size optimized field intensity with 
yellow showing super-oscillating regions. 



local vector (gradient of the phase) larger than the spectral bandwidth of the original field. 
To visualize this effect, in the case of OEi spot size optimized beams, we have calculated the 
spectral density of the radial wave- vector for the smallest planar spot 1 16]. As shown in Fig. [3] 
this spectral density clearly identifies a spectral bandwidth (white background in Fig. [5]). Re- 
gions of the beam which exhibit locally larger wave- vectors than the ones supported by this 
spectral band width correspond to super-oscillating regions. The local wave vector is defined 
as d r arg(w(r)) where arg(w) defines the phase of the analytical signal u. In this particular case, 
we observe that super-oscillations occur in the dark region of the beam. Additionally, when the 
ROI is large compared to the Gaussian beam waist wo, there are no super-oscillating regions. 
These only appear when the beam starts to be squeezed. 

2.2. Smallest focal spot using Bessel beams 

The paraxial approximation employed above in the case of LG beams can be used to describe 
sub-diffracting beams but breaks down when beams are tightly focused. As a consequence we 
must consider full vectorial solutions of Maxwell's equations. Here, we have chosen Bessel 
beams as a basis-set and determined the superposition of Bessel beams which minimized the 
spot size in a planar finite ROI. Note that the problem of the finite intensity of Bessel beams ifTTl 
is easily circumvented here due to the finite ROI size considered. The monochromatic electric 
vector field of the vectorial Bessel beam may explicitly be expressed as fT8) 

E = Eoexp(i£(j) + ik t z)) ^(ae x + Pe y )Ji(k t r) 

+ lr ^ a + ^ ex PH0)^-i (V) - (a - *P) exp(z0)/ m (k t r))e z j (8) 

where k t = ko sin(0) and k z = &ocos(0) are the transversal and longitudinal wave vectors with 
the characteristic cone angle of the Bessel beam. e x , and e z are the unit vectors in the 
Cartesian coordinate system. The parameter £ corresponds to the azimuthal topological charge 
of the beam while a and f3 are associated with the polarization state of the beam. The magnetic 



field H was deduced according to Maxwell's equations. Figure]?] shows a comparison between 
the Airy disk, the Bessel beam and the OEi optimized spot considering a numerical aperture 
of NA= 0.1. As in the case of the LG beams, squeezing the focal spot is accompanied by side 
bands and a loss in efficiency shown by the Strehl ratio (see Fig. [5}. 




Fig. 4. Intensity cross sections: (a) Airy disk for the maximum numerical aperture consid- 
ered NA= sm(6 max ) =0.1. The yellow dashed circle shows the position of the smallest 
zero-intensity circle taken as the ROI inside which the spot size is calculated. The spot size 
is normalized to the spot size of the reference Bessel beam, (b) Reference Bessel beam 
corresponding to the largest cone angle 6 max . The spot size of the reference Bessel beam 
is denoted as wb- (c) OEi spot size optimized beam for a superposition of Bessel beams 
(0 G [0, Qmax]) for a large ROI highlighted by the dashed yellow circle. Strehl ratio: 2%. (d) 
OEi spot size optimized beam for a small ROI. Strehl ratio: 0.2%. The gray-scaled region 
shows the sidebands while the color range the ROI. Notice that the two scales are different. 




Fig. 5. (a) Relative spot size Ar/w^ of the Bessel beam superposition as a function of the 
relative ROI radius R/Rb- The spot size wb and the ROI radius Rb are associated with the 
reference Bessel beam shown in Fig.Qb), where the ROI is indicated as dashed circle. For 
comparison, the red dot indicates the location of the reference beam in the Ar/wB vs. R/Rb 
plot, (b) Strehl ratio vs relative ROI radius R/Rb- 

3. Experimental OEi 

3.1. Experimental implementation of the OEi concept 

To perform an experimental OEi optimization we have used the following setup: A HeNe laser 
beam is expanded and subsequently amplitude modulated by a spatial light modulator (SLM) 
display operating in conjunction with a pair of crossed polarizers. Analogously to a liquid crys- 
tal display on a computer or laptop monitor, the liquid crystal SLM display rotates the polariza- 
tion of the incident light by an angle depending upon the voltage applied to the display pixels. 
The amplitude modulated beam is then imaged onto a second SLM display through a pair of 
lenses. This second SLM display along with a subsequent Fourier lens and aperture served to 
modulate the phase of the laser beam in the standard first order configuration [20]. The field 
modulations of interest were encoded as RGB images where the blue channel represented the 



amplitude and the green channel the phase modulation. The SLM controller extracted these 
information and applied the two channels to the respective panel. We have performed calibra- 
tion measurements to ensure that both the amplitude and phase modulation exhibited a linear 
dependence on the applied 8-bit color value between and 255. A CCD camera allowed us to 
record images of laser fields in the Fourier plane of lens 5. 
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Fig. 6. Experimental setup. FP = focal plane, L = Lens. Focal widths: f\ = 50 mm, 
f2 = 500 mm, fo = j\ = 400 mm, f$ = 1 m. Laser: JDS Uniphase HeNe laser, P max = 
10 mW, A = 633 nm, SLM: Holoeye HEO 1080 P dual display system, resolution = 
1920 pixel x 1080 pixel, display size = 1 in x 0.7 in. CCD camera: Basler pilot piA640- 
210gm, resolution = 648 pixel x 488 pixel, pixel size = 7.4 um x 7.4 um. 



The experiment consist in determining the OEi in the CCD camera plane whilst shaping and 
superimposing the test fields Ej in the SLM planes. In the following, we indicate the plane 
of interest by a z-coordinate along the optical axis where z = z\ and z = zi refer to the SLM 
and CCD camera plane, respectively. According to this convention we shape a set of test fields 
Ej{zi) = Aj(zi)e l fy( Zl ' both in amplitude Aj and phase <j>j in the SLM plane, and the associated 
intensities Ij{zi) \Ej(z2)\ 2 are detected in the CCD camera plane. The amplitudes Aj(zi) 
were determined from these intensities by simply taking the square root. We used the three- 
step phase retrieval algorithm described in Ref. I2H to retrieve the phase modulations tyj{zi). 
The determination of the phase and amplitude of the beam in the CCD plane allows us to 
numerically vary the ROI without redisplaying the test fields. Using these fields, the IO and 
S SO are determined according to Eqs. ([5j and ([6]), respectively. 

During the course of our experiments we verified the linearity of our optical system by per- 
forming a comparison between what we term the "experimental superposition (Exp-S)" and 
the "numerical superposition (Num-S)". The Exp-S refers to the case where the set of OEi op- 
timized superposition coefficients ai is used to encode the optimized superimposed field onto 
the SLM. The CCD camera then detected the intensity /Exp-sfe) corresponding to this en- 
coded optimized field. The Num-S utilizes the fields Ej{z%), which were individually measured 
to assemble the OEi operators, in order to numerically determine the intensity distribution as 

^Num-sfe) <* \H?=iajEj(z2)\ 2 . Crucially, linearity is verified if 7 E x P -sfe) = ^Num-sfe). This is 
indeed observed in our experiments as demonstrated in the following subsection which features 
a comparison of experimental and numerical intensity distributions. 



3.2. Results and discussion 



In our experiments, we used N = 11 non overlapping amplitude ring masks with a constant 
phase modulation as fields of interest Ei{z\). After propagation through the Fourier lens 5 (see 
Fig. [6]) the resulting fields Ei{i2) form a set of Bessel beams. Figure [TJa) shows the largest 
ring modulation encoded onto the SLM with the resulting Bessel beam shown in Fig. (b). As 
this particular Bessel beam comes along with the highest NA compared to the Bessel beams 
created with smaller ring modulations, the beam shown in Fig. (b) exhibits the smallest central 
spot of all beams realized in our experiments. The spot size of the Bessel beam featuring the 
smallest core is denoted as wb and used as reference for the measurements presented below. 
For comparison Figure [TJc) depicts a circular aperture which is encoded onto the SLM in order 
to observe the Airy disk (see Fig. (d)). The spot size of the Airy disk is approximately 1 .5 times 
larger than the core of the reference Bessel beam as expected l22l . 
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Fig. 7. SLM encoded field modulations and resulting beam profiles, (a) Ring mask RGB 
image as encoded onto the dual panel SLM. (b) Associated Bessel beam created in the 
CCD camera plane, (c) Aperture RGB image as encoded onto the dual panel SLM. (d) 
Associated Airy disk as detected by the CCD camera. The yellow bar in (b) represents 2 
times the spot size w# of the Bessel beam's central core, w in (d) is the spot size of the Airy 
disk. 



The results of the performed OEi spot size minimization are shown in Fig. [8] for different 
sizes of the ROI. To begin with, the comparison of the Num-S intensity distribution /Num-sfe) 
(top row) and the Exp-S intensity distributions ^Exp-sfe) (bottom row) clearly reveals good 
agreement and thus verifies the linearity of our optical system as elucidated above. For com- 
pleteness, the central row shows the Exp-S superposition in RGB format as encoded onto the 
SLM. The color code features a blue channel representing the amplitude modulation from 
(black) to 1 (blue) and a green channel corresponding to phase modulations from (black) to 
2tc (green). Next, we conclude from the measured relative spot size w/wb that the spot size 
decreases if the ROI size is reduced. The reduced spot size is achieved at the expense of the 
spot intensity which is redistributed to a ring outside of the ROI similar to the theoretical results 
presented in section 22_ and Fig. [4] Referring to the Exp-S data, for R = 7 pixel the spot size 
is reduced to 72% of the size of the reference Bessel beam's core and even further to 50% for 
R = 4 pixel. The latter result is somewhat vague, though, due to the low spot intensity which 
may be truncated by the sensitivity threshold of the CCD detector and thus may appear smaller. 
However, our experimental results overall clearly verify the OEi concept applied to spot size 
minimization. Moreover, the results strongly suggest that the OEi optimization may indeed 
squeeze spots to the subdiffractive regime since the optimal superposition of Bessel beams not 
only beats the Airy disk but also the reference Bessel beam diameter. 



4. Applicability of the OEi method to scattering interactions 

In this section, we demonstrate, on the basis of a numerical study, how the OEi method can 
equally be applied in some scattering interaction processes. Indeed, the OEi method presented 
above is applicable to free space propagation and can be directly extended to linear scattering 
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Fig. 8. Experimental OEi spot size minimization. Top row. /Num-sfe) f° r different ROI 
radii in pixel as indicated in the top left corner of all graphs shown. The ROI is exemplary 
indicated as a dashed ring in the left hand side intensity distribution. The number in the 
bottom left corner represents the spot size w in units of the reference spot size w# . Central 
row. Optimized experimental distribution as RGB encoded onto the SLM. Bottom row. 
Intensity distributions /Exp-sfe)- The relative spot size w/w# is indicated in the lower left 
corner. 

processes where the optical interaction is expressed as a quadratic form of the field. As in the 
case of the smallest spot operator, the OEi of these light matter interactions can be used to deter- 
mine the electromagnetic field profile delivering the largest or the smallest interaction strength. 
In this section, we show two numerical examples illustrating the cases presented in Table [T] 
The numerical modelling is performed using a finite element method (Comsol) and solving the 
fully vectorial monochromatic Maxwell's equations in 3D. The structures considered here are 
embedded in a larger computational domain surrounded by perfectly matched layers. Figure [9] 
shows the electric field amplitude |E| for the different cases considered. To implement the OEi 
method, we determine the matrix operator with the help of equation ([2]). Here, we use angular 
spectral decomposition |fT9l of the incident light field corresponding to a numerical aperture of 
NA=0.8. 

In the first example (Fig.[9]a-c), we use the intensity operator associated to the measure 
to determine the largest transmission through a sub- wavelength aperture (diameter=200nm) in 
a thin layer of silver (thickness=200nm). The incident light field considered is linearly polarised 
and the transmission is determined across the output surface of the aperture. 

The electric amplitude |E| of most efficient transmission OEi is shown in Fig. [9^ illustrating 
scattering of the aperture. The OEi on its own in Fig.[5j). The transmission enhancement factor, 
with respect to the tightest Bessel beam achievable for a numerical aperture of NA=0.8 (Fig. 
[9};), is 2.1 and 1.55 with respect to the Airy diffraction limited disk with the same numerical 
aperture. 

In the second example, we use the optical momentum operator associated with the quadratic 
measure m( F ' Uz ) as defined in Tablehl The momentum OEi with the largest positive eigenvalue 
corresponds to the field profile (FigH9]e-f) giving the largest optical force on the microparticle. 




Fig. 9. (a-c) Cross section plot of the electric field amplitude, |E|, for a sub- wavelength 
aperture (diameter=200nm) in a thin layer of silver (thickness=200nm, refractive index 
n = 0.12 — 3.7/ at a wavelength A = 600nm) illuminated from below. The yellow lines 
represent the boundary of the structure, (a) Intensity OEi ensuring the largest transmission 
(transmission enhancement factor 2.1 with respect to the tightest Bessel beam and 1.55 
with respect to the Airy disk illumination), (b) Incident intensity OEi without the structure, 
(c) Tightest Bessel beam illumination, (d-f) Electric field amplitude, |E|, in a cross section 
for a high refractive index (n = 1.8) microparticle (diameter=800nm) illuminated from 
below with a wavelength (A = 504nm). (d) Momentum OEi ensuring the largest momentum 
transfer (enhancement factor 49.3 with respect to the plane wave and 1.33 with respect the 
Airy disk illumination), (e) Incident momentum OEi without the structure, (f) Plane wave 
illumination. 

Figure[9|l shows the field amplitude |E| of this OEi on its own and scattering of the microparticle 
(Fig.[9E). The optical force enhancement factor, with respect to the plane wave illumination 
(Fig.[9p, is of 49.3 and of 1.33 with respect to the Airy diffraction limited disk with the same 
numerical aperture. 

5. Discussion and Conclusion 

We have experimentally and theoretically demonstrated an approach based on optical eigen- 
modes that enables the minimization of the free space spot size of a beam. Using full vectorial 
simulations in 3D, we have shown how the OEi approach can be used to optimise light-matter 
scattering interactions in the case of transmission through sub-wavelength apertures and optical 
forces on micro-particles. The generic nature of our approach means that it can be applied to 
other cases where the measure has a quadratic form and the propagation is linear. In the present 
paper we have verified the rigor of the method by demonstrating the experimental spot size 
operator and intensity operator optimization using Laguerre-Gaussian and Bessel light modes 
using a dual SLM to implement the technique. Future work will aim to extend this method to 
optimise the size and contrast of optical dark vortices, the linear Raman scattering or the fluo- 
rescence of different samples, the optical trapping force, and the angular momentum transfer in 
optical manipulation. 
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Annex: Dimensionality study 

In this appendix, we verify the convergence properties of our optical eigenmode method as a 
function of the dimensionality of the probing base used. Here, we consider a plane wave basis 
corresponding to the angular spectral decomposition of any incident field. The propagation is 
modelled using direct Fourier transform optics between the SLM plane and the target plane. 
Additionally the phase front of the incident beam is randomly changed to simulate the effect of 
high aberrations within the optical train. The dimensionality for the optical eigenmode method 
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Fig. 10. (a) Comparison between the enhancement factor achieved using standard phase 
front correction techniques (in black [23 24 1), the phase front correction from the intensity 
eigenmode (in red) and from the smallest spot size eigenmode (in blue), (b) Normalised 
beam spot size for the standard phase front correction techniques (in black [23 ( 24 1) and 
from the smallest spot size eigenmode (in blue). COq corresponds to the spot size of the Airy 
disk. 



corresponds to the number of plane waves, N, taken into account, whilst retaining a fixed nu- 
merical aperture. The convergence behaviour is compared to standard phase front correction 
methods which can also be used to achieve focalised spots in the case of highly aberrated light 
fields. More precisely, these methods are based on the variable partitioning of the SLM to cre- 
ate N beamlets whose phases are individually changed such that all beamlets constructively 
interfere in the focal target point l23ll24l . This approach delivers a final correction phase mask 
that is able to correct for aberrations in the SLM incident light field additionally to the pre- 
correction of the propagation aberration between the SLM and the target focal plane [24]. This 
correction mask delivers an enhanced focal intensity rj shown in black in Fig. [TOfr for both 
standard methods E3 For comparison purposes, we define a phase front correction mask 
using the phase part of the intensity operator eigenmode (red in Fig.flOh) which shows that, for 
all three methods, the enhancement factor scales linearly with the dimension N. 

Further, Fig 10 3 shows that the standard approaches cannot beat the diffraction limit, given 
by the Airy disk. In stark contrast, the smallest spot size eigenmode delivers sub-diffraction 
over almost the whole dimensional range considered with increasing efficiency (blue curve in 
Fig.flOk) for increasing dimensionality. 



